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A method for determining the positions of hologram dislocations relative to the optical axes 
of entangled Laguerre Gaussian modes is proposed. In our method, the coincidence count rate 
distribution was obtained by scanning the position of one of the holograms in two dimensions. 
Then, the relative position of the hologram dislocation was determined quantitatively from the 
positions of the minimum and maximum coincidence count rates. The validity of the method was 
experimentally verified, and in addition, an experiment demonstrating the violation of the Clauser- 
Horne-Shimony-Holt inequality was performed using the well-identified optical axes of the entangled 
modes. 
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I. INTRODUCTION 

The quantum states of photons in Laguerre-Gaussian 
(LG) modes have been attracting a lot of attention re- 
cently. It is well known that photons in LG modes have 
orbital angular momentum, just as circularly polarized 
photons have spin angular momentum. However, there 
is a big difference between these two physical properties. 
In principle, aniVxiV Hilbert space with arbitrary large 
N can be created from a single photon in an LG mode, 
which is much larger than the 2x2 Hilbert space ob- 
tainable with polarization. Such quantum states, called 
quNits-can be used to improve quantum information pro- 
tocols PJ, Ll L| • It is a l so interesting to study the nature 
of entanglement in larger Hilbert spaces using single pho- 
tons in LG modes. 

In a pioneering paper 4], Vaziri et. al. confirmed 
experimentally that the angular momentum of photon 
pairs are conserved through the parametric down conver- 
sion process and that the generated pairs can be repre- 
sented as a superposition of states in different LG modes. 
As a logical conclusion, they claimed that the gener- 
ated pairs are entangled in LG modes. Later, the same 
authors reported a violation of Clauser-Horne-Shimony- 
Holt (CHSH) inequality in three dimensions 5]. Re- 
cently, Langford et. al. performed quantum state tomog- 
raphy with their entangled photonic qutrits and qubits 
and confirmed entanglement quantitatively 6]. Concen- 
tration of entanglement in LG modes 7], and produc- 
tion, transmission, and reconstruction of qutrits |8| have 
also been demonstrated using entangled photons in LG 
modes. 

In those experiments, holograms with a dislocation in 
the center were used as LG mode converters 0. When 
the optical axis of the light beam is 'exactly' located on 
the dislocation, a state in a LG mode is converted to one 
in a Gaussian mode, and is detected by a photon counter 
through a single mode optical fiber. In contrast, a su- 
perposition state of different LG modes is converted to 



a Gaussian mode when the optical axis is located at an 
appropriate position relative to the dislocation. Thus, it 
is crucial to be able to determine the position of the dis- 
location relative to the optical axis. However, it is not a 
trivial matter for the case of entangled LG modes, typi- 
cally generated via spontaneous parametric down conver- 
sion (SPDC), because the optical axis of the signal beam 
have to correspond with the optical axis of the idler beam 
in accordance with the phase matching condition. 

In this paper, we present a method for determining 
the position of the hologram dislocation relative to the 
optical axis of entangled LG modes. In our method, we 
obtain the coincidence count rate distribution by scan- 
ning the position of the hologram relative to the idler 
mode position in two dimensions. The relative position of 
the hologram dislocation in the signal beam can then be 
determined quantitatively from the location of the mini- 
mum and maximum points of the coincidence counts [l^. 
We also found that the obtained two-dimensional maps 
of the coincidence count rates exhibited intuitively a dis- 
tinctive feature of entanglement, namely, correlations in 
different non-orthogonal basis sets. From the experi- 
mentally obtained two coincidence distributions, we suc- 
ceeded in showing, using the well-defined optical axes, 
that violation of the CHSH inequality is a signature of 
entanglement in LG modes. 

This paper is organized as follows. In section II, we 
explain how LG modes are converted by reflection-type 
hologram gratings. In section III, we first explain how we 
can identify the measurement basis states by considering 
the time-reversal process, and then we describe the pro- 
posed method for determining the position of the optical 
axis relative to the hologram dislocation. In section IV, 
we describe our experimental setup. In section V, the 
procedure used for determining the relative position of 
the hologram dislocation to the optical axis from exper- 
imental measurements is described in detail. In section 
VI, we report the violation of the CHSH inequality, us- 
ing the well-defined optical axes of entangled LG modes. 
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Then, we conclude this paper in section VII. 

II. CONVERTING LG MODES WITH 
HOLOGRAMS 

When a pump light in LGo mode is used, the quantum 
state of photon pairs generated through the process of 
SPDC is written as follows [Toj ] . 

+ CO 

|$) = a m |-ra) A |ra)B- (1) 

171— — OO 

Where \m) is a one-photon state of a LG m mode, and 
a m is the probability amplitude. Here we neglected the 
radial index of the LG modes since the radial index was 
not distinguished in our experiment. 

In order to detect photons in a specific superposition 
of states of LG modes, an optical setup with a hologram 
and a single mode fiber (SMF) 11] shown in Figda) was 
used. Here we assume that a reflection type hologram 
is used and that it has a one-pronged dislocation at its 
center We also adopt the definition of diffraction 

order number given in Fig. IHb). When the central axis 
of the incident beam in a LG n mode is centered on the 
dislocation of the hologram, the m-th order diffracted 
mode is converted to a LG_ n _ m mode. Thus, a photon 
in the LGi mode in the incident beam will be converted 
to a LGo mode in the —1st order diffracted beam. In 
other words, the photons which are not in the LGi mode 
will not be converted to LGo mode. The SMF, which is 
positioned in the —1st order diffracted beam, functions 
as a filter that blocks all the LG modes except for the 
LGo mode. Thus, by considering the time-reversal of 
this process, we can conclude that the photons which are 
detected after the SMF are those which were in the LGi 
mode in the incident beam. 

In order to measure the photons in a superposition 
state of different LG modes, the position of the dislo- 
cation needs to be shifted from the central axis of the 
incident mode. We define r as the distance between the 
dislocation in the hologram and the central axis of the 
optical beam, and as the azimuthal angle illustrated 
in Fig. [21 Suppose that photons in a Gaussian mode 
are incident onto the hologram. Since the hologram is 
designed so that the position of the phase singularity of 
the diffracted beam is located on the dislocation in the 
hologram, the diffracted photon will have a superposi- 
tion of states in the LGo and LGi modes, given, to a 
good approximation, by the following expression |l7j . 

, where uo is the beam radius. Note that the diffracted 
photons will also be in other modes (i.e. LG_i, LG2, 




FIG. 1: The experimental set up for detecting photons in 
LGi mode. If the dislocation of the hologram is located on 
the central axis of the incident beam of LGi mode, then the 
-1st order diffracted beam will be in LGo mode, which couples 
into the SMF. 

(a) Position of hologram 

r \ dislocation 




Position of optical axis 

(b) 




FIG. 2: An experimental set up for detecting photons in a 
superposition of LGo and LGi modes. When the dislocation 
of the hologram is located at the position of the phase sin- 
gularity of mode, the -1st order diffracted beam becomes the 
LGo mode which couples into the SMF. 



and so on), however, here we neglect those terms since 
the amplitudes of those states are small [HI . 

When an incident photon is in a pure state of one or 
two LG modes, it is easy to determine the position of the 
dislocation in the hologram relative to the optical axis 
of the modes. For instance, the intensity distribution 
of the diffracted beam directly gives the position of the 
phase singularity point, which is where the intensity of 
the beam is zero. However, this is not the case for en- 
tangled photon pairs in LG modes. This is because the 
quantum state of one photon in the pair is not a pure 
state, but a mixed state of the various LG modes, and 
the intensity distribution (the spatial distribution of the 
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SMF detector 

FIG. 3: An experimental set up for the measurement of cor- 
relations between LG modes in path A and path B. 

single count rates) does not show any singularity. 

In the next section, we propose a method for estimat- 
ing the position of the hologram dislocation relative to 
the optical axis of the entangled photon pair, which in- 
volves taking coincidence count probabilities by scanning 
the hologram position in two dimensions. 



III. DETERMINING THE OPTICAL AXIS 
POSITION USING COINCIDENCE COUNT 
RATES - THEORY 

In the following calculation, we assume that the signal 
photons and idler photons propagate along the optical 
paths A and B, and then the -1st order and 1st order 
diffraction beams are coupled to two SMFs, respectively 
(Fig. 03). The position of the dislocation relative to the 
optical axis for each path are (va,0a) an d (re, 0b). 

In order to calculate the coincidence count probabili- 
ties, we first need to know the measurement basis states 
for the given optical set up. In other words, we have to 
identify the 'state', which is converted into a Gaussian 
mode as a result of diffraction by the hologram. This 
can be found by considering the time-reversal process 
discussed in the previous section. 

We have already stated that a Gaussian state is con- 
verted into the state given by eq.(2). Noting that 
exp(— iff)\m) is converted into exp(i0)(— l) m |ra) by the 
time-inversion operation, we find that the superposition 
state given by, 

I = «- tf ^^|0> + ^^H>. (3) 

is converted to the Gaussian state |0) in the -1st order 
diffraction beam. Therefore, eq.© is a measurement ba- 
sis state for path A in our experimental setup. 



Using a similar argument, we also find that 

is converted to the Gaussian state |0) in the 1st order 
diffraction beam. Similarly, eq.Q is a measurement basis 
state for path B in our experimental setup. 

From eq. {IJ , (|3J) , and (@J , the coincidence counts prob- 
ability P for the entangled photon pair is given by 

= \A(-(r A ,e A )\B(Hr B ,e B )\<s>)\ 2 

■ |2 ^A r B + ^r A r B oo 2 cos(-0 A + Ob + 6) + co 4 
W (2r\ + ^ 2 )(2r|+^ 2 ) [) 

, where we assume that — a and ao = ae %5 . 

This equation suggests that, as a result of entangle- 
ment, the coincidence probability has a maximum and 
a minimum (zero) at particular relative positions of the 
hologram in path A, for a given relative position (r#, Ob) 
of the hologram in path B. Actually, P has a maximum 
value \a\ 2 when (ta^Oa) — ( r B,5 + 0b), and a minimum 
value when (ta, 0a) = (^ 2 /(2re), tt+S+Ob)- Such min- 
imum and maximum positions can be found by taking co- 
incidence count probabilities while scanning the position 
of hologram A. Note that the maximum and minimum 
positions are always opposite to each other relative to the 
position of optical axis (ta — 0) between them. Conse- 
quently, the distance d between the maximum position 
and the minimum position is given by 



In order to check the validity of this method, we per- 
formed the experiment described in the next section. 

IV. EXPERIMENTAL SETUP 

The experimental set up used is shown in Fig. 21 We 
used a /3-barium borate (BBO) crystal cut for Type I 
phase matching condition having dimensions 5 (height) 
x 5 (width) x 3 (thickness) mm. It was pumped by an 
argon- ion laser, having a wavelength of 351 nm, a beam 
radius of 1.2 mm and a power of 150 mW. The generated 
signal and idler photons were focused by a plano-convex 
lens (/ = 300 nm) and the hologram was positioned in 
the focal plane. The hologram plate was made by fabri- 
cating the structure on thin polymer layer on glass sub- 
strate using an electron beam writer, and then coating 
with gold. The hologram was 2 mm in diameter. The 
-1-th (1-th) order diffraction beam in path A (B) was 
collimated by a plano-convex lens (/ = 200 nm) and 
filtered by a narrow band-pass filter (IF) having a cen- 
ter wavelength of 702 nm and a FWHM bandwidth of 
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FIG. 4: Experimental setup used to determine the hologram 
position 



Measured coincidences Simulation 




Horizontal position [urn] Horizontal position [lliii] 

FIG. 5: A comparison of measurement data with calculation 
data 

4 nm. The photons were then coupled to SMFs via ob- 
jective lenses, and counted by single photon detectors 
(AQR-FC, Perkin Elmer). The single counting rates and 
coincidence counting rates were measured by a photon 
counter (SR-400, Stanford Research System). 

The single counting rates for both paths A and B were 
maximized by adjusting the positions of the respective 
coupling lenses. While keeping the position of the holo- 
gram B fixed, we measured the coincidence counting rates 
while scanning in two dimensions the position of the holo- 
gram A, which was positioned in the focal plane of the 
first convex lens. Typical single counting and coincidence 
counting rates were 5000 [cps] and 100 [cps], respectively. 
The vertical and horizontal positions were shifted by a 
step of 150 /im for a grid made up of 14 x 14 points. It 
took 40 minutes to obtain a complete map of the coinci- 
dence count rates. 



V. DETERMINING THE OPTICAL AXIS 
POSITION USING COINCIDENCE COUNT 
RATES - EXPERIMENT 

Fig. El (a) is a contour plot of the coincidence count 
rates obtained by scanning the position of hologram A for 
a certain position of hologram B. The horizontal and ver- 



tical axes represent the horizontal and vertical position 
of hologram A relative to the optical axis, respectively. 
The white area indicates the location of the hologram 
where the coincidence rates were highest. Note that the 
measured coincidence counting rates were affected by the 
position-dependent diffraction efficiency of the hologram 
we used; the measured efficiency of hologram A decreased 
as the distance from the dislocation increased. The coin- 
cidence rate was almost zero when ta is large ( > 1mm 
). This may be because the incident beam was outside 
of the hologram area (2 mm in diameter). 

The measured coincidence count distribution had a 
maximum at (x,y) = (—50/im, 200/im) and a minimum 
at (x,y) = (—50/im, —400/im). As discussed in the pre- 
vious section, the position of the optical axis in path 
A should be located half-way between these two points. 
Therefore, b + S is estimated to be tt/2. 

In order to estimate the other parameters, tb and uo, 
we tried to simulate the coincidence probabilities using 
eq.JJJJ). The position-dependent diffraction efficiency of 
the hologram, which was experimentally measured using 
a HeNe laser beam (632 nm) while scanning the focus 
position, was multiplied by P(ta^a) in eq.JSJ). In Fig. 
Eta)? the distance between the maximum point and the 
minimum point was 600 /im. Thus from eq. ©, tb + 
uo 2 /(2vb) = 600/im. We performed several calculations 
for different uo and tb which satisfy this relation and 
found that the distribution of the high coincidence region 
is well approximated by uo = 400/im and r# = 200/im. 
The result of simulation is shown in Fig. Et D ). The 
estimated uo is consistent with our optical setup. 

In order to determine the remaining parameters, we 
shifted the position of hologram B and measured the 
maps of coincidence count rates by scanning hologram 
A. The results for nine different positions of hologram B 
are shown in Fig. El From all of these maps, Fig. Et e )> 
which was measured by shifting hologram B by 200 jam 
upwards, exhibits a particularly symmetric pattern. For 
this condition, the minimum point (uo 2 / (2r b) tt + S + Ob) 
should be far from the center and thus tb should be close 
to zero. That is, the dislocation of hologram B is almost 
on the optical axis. As a result, for the original hologram 
B position in Fig. El fB was confirmed to be 200 /im and 
0b was determined to be —tt/2. We have already seen 
that 0b -\-S = tt/2. Consequently, S is calculated to be n. 

The simulated coincidence count probabilities calcu- 
lated following the same procedure as that of Fig. Etb) 
are shown in Fig. {7\ In the simulation, the maximum 
point appears on the opposite side of the origin from the 
minimum point, and both points move as the relative 
position of the hologram B changes. This is due to the 
non-classical correlations in the different non-orthogonal 
basis sets. These characteristics are observed in the ex- 
perimental results (Fig. El- 
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FIG. 6: Maps of the coincidence count rates obtained by scan- 
ning hologram A with nine different positions of hologram B. 
In (e) , the dislocation of hologram B is almost on the optical 
axis. The reference position of hologram B was taken to be 
(tb,0b) = (0/ira, 0). The positions of hologram B relative 
to this position are (a): (200/im,37r/4), (b): (200/xm,7r/2), 
(c): (200/xm,7r/4), (d): (200/xm,7r), (f): (200/xm,0), (g): 
(200/xm,7r/2), (h): (200/ira, -tt/2), (i): (200/xra, -7r/4). 
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FIG. 7: The simulated coincidence count probabilities for nine 
different positions of hologram B. The positions of hologram 
B (r B ,0 B ) are (a): (200/xm, 3tt/4), (b): (200/im,7r/2), (c): 
(200/xm,7r/4), (d): (200/xm,7r), (e): (0/xm,0), (f): (200/xm,0), 
(g): (200/xm,7r/2), (h): (200//m,-7r/2), (i): (200/zm ,-tt/4). 



VI. VERIFICATION OF ENTANGLEMENT 

We also tried to verify the non-local correlation us- 
ing the CHSH inequality^!, using the well-defined op- 
tical axes of the entangled modes. In the original 
discussion [l^. Clauser and co-authors assumed that two 
photons are sent along two opposite paths A and B, and 
the polarizations of those photons were analyzed using 
polarization analyzer in both paths. Those polarization 
analyzers were characterized by a parameter and they 
discriminate the photons into two orthogonal channels, 
channel and 0, where the photons are detected. Under 
the 'fair sampling assumption' 0, a correlation function 
E(6a,0b) is defined as follows: 



E{0a,0 b ) 

C(0 a ,0 b ) + C^0b~) 



C(0 a ,0b)-C(0 a ,0b) 



C(0 A , b ) + C(6 A , B ) + C(0 A , 0b) + C(6 A , B ) 

(7) 

Here, C(Q A , 0b) is the coincident count rate of photons 
in channel A and B • Then, the CHSH inequality, 

\S\ = \E(0 A ,8 B )-E(8' A ,9 B )+E(8 A ,8' B )+E(0' A ,8' B )\ < 2, 

(8) 

is valid for any local hidden variable theory. 

In order to verify the violation of the CHSH inequality, 
we measured coincidence counts by changing the relative 
azimuthal angle Ai 0b of the holograms A and B, while 
keeping the relative distances constant, r A — tb = r. 
In this case, the coincidence probabilities are given as a 
function of A and ^ as follows, 



P(0 A ,0 B ) = a [ 



, 4r 4 + u 4 + 4r 2 u 2 cos(5 + A - B ) 
' (2r 2 +cj 2 ) 2 



. (9) 



In our experiment, we were just concerned with the two 
dimensional Hilbert space of each path A and B, since 
our experimental setup just measured the components of 
|0) and |1) in path A, and of |0) and | — 1) in path B. 
In this case, we can substitute photon detection of the 
parameter 1 - = + ir for that in the channel 0. That is, 
the coincidence rate C(0 Al 0b) in eq. (JSJ) can be replaced 
byC(0i,0 B )0 . 

E(9 A ,9 B ) 

_ C(8 A ,0 B ) + C(8i,8^) - C(8i,8 B ) - C(8 A ,8^) 



C(0 A , B ) + C{8\, 0%) + C(0Ji,0B) + C(6 A ,6 



(10) 

If it is assumed that the channel loss is independent of the 
angle then we can calculate the correlation function 
E(0 Ai 0b) by replacing the coincidence rates C(0 a ,0b) 
by P(0 Ai 0b) in eq. (fTU|). When we set the parameters 
A = -tt/4, 0' a = tt/4, b = -tt/2, 0' b = and r = 
uj/V2, we obtain \S\ = 2^/2 ~ 2.828(> 2) 
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TABLE I: Coincidence count rates measured over 5 seconds. 



We performed experiments with the parameters r = 
200/im, A = 0' a = tt/4, b = -tt/2, 0' b = 0. The 

obtained coincidence rates over 5 s accumulation period 
are shown in Table 1. Note that these data are measured 
independently from the data shown in Fig. The ob- 
tained |5| was 2.127, which is larger than 2, indicating 
that the CHSH inequality has been violated |2Cj . 

The reason why | aS' | was smaller than the maximum 
value 2.828 seems to be mainly because the value for 
r we used, r — 200/im, was smaller than the optimum 
position of r = 283/im. In fact, 1*51 is calculated to be 
2.26 for r = 200/im, which is close to the experimental 
result. 



VII. CONCLUSION 

In conclusion, we have proposed a method for deter- 
mining the position of the hologram dislocations rela- 
tive to the optical axes of entangled Laguerre Gaussian 
modes. In our method, the coincidence count rate distri- 
bution is obtained by scanning the position of one of the 
holograms, and then the relative position of the hologram 
dislocation is identified quantitatively from the positions 
of the minimum and maximum coincidence count rates. 
We tested this method experimentally with the entan- 
gled photons generated thorough the SPDC process and 
succeeded in determining the positions of the hologram 
dislocations relative to the optical axes of each mode. Fi- 
nally, an experiment investigating the CHSH inequality 
was performed using the well-determined optical axes of 
the entangled modes, and the obtained value \S\ = 2.121, 
which is slightly smaller than the theoretically predicted 
value of 2.26, demonstrates violation of the CHSH in- 
equality. 

In the theory part, we neglected the radial index of 
the LG modes and derived eq. (5). It should be noted 
that the proposed method based on eq. (5) is useful as is 
shown in the experimental demonstration. However, the 
further studies on the effect of radial indexes and other 
'outer modes |g' on the entanglement in LG modes as 
well as new experimental method to avoid those effects 
seem interesting and important. 

The authors would like to thank J. Arlt for provid- 
ing us the reference H. F. Hofmann, K. Tsujino, R. 
Okamoto for their discussion, M. Moriya for gold coating, 
H. Ohminato, T. Yonemura for their technical support. 
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